A STABILITY INDEX FOR DETONATION WAVES 
IN MAJDA'S MODEL FOR REACTING FLOW 



Gregory Lyng and Kevin Zumbrun 
October 6, 2000; Revised: May 28, 2003 

Abstract. Using Evans function techniques, we develop a stability index for weak 
and strong detonation waves analogous to that developed for shock waves in [GZ,BSZ] , 
yielding useful necessary conditions for stability. Here, we carry out the analysis in 
the context of the Majda model, a simplified model for reacting flow; the method is 
extended to the full Navier— Stokes equations of reacting flow in [Ly,LyZ]. The result- 
ing stability condition is satisfied for all nondegenerate, i.e., spatially exponentially 
decaying, weak and strong detonations of the Majda model in agreement with nu- 
merical experiments of [CMR] and analytical results of [Sz,LY] for a related model 
of Majda and Rosales. We discuss also the role in the ZND limit of degenerate, sub- 
algebraically decaying weak detonation and (for a modified, "bump-type" ignition 
function) deflagration profiles, as discussed in [GS.1-2] for the full equations. 



Section 1. Introduction 

In one-dimensional, Lagrangian coordinates, the Navier-Stokes equations of re- 
acting flow for a one-step reaction may be written in the abstract form 



u t + f(u) x = (B(u)u x ) x + kq(p(u)z, 
z t = (D(u, z)z x ) x - kcp(u)z, 



where u, f, q G K n , B G R nXn , z, k, D, ip G R 1 , and k > (model (8.81) of 
[Z.3] with particle velocity set to zero). Here, vector u comprises the gas-dynamical 
variables of specific volume, particle velocity, and total energy, and z measures mass 
fraction of unburned reactant: more generally, "progress" of a single reaction in- 
volving multiple reactants. The first equation thus models kinematic and the second 
equation reaction effects. The function <p(u) is an "ignition function", monotone 
increasing in temperature, and usually assumed for fixed density to be zero below a 
certain ignition temperature and positive above. The vector q comprises quantities 
produced in reaction, in particular heat released. The coefficient k corresponds to 
reaction rate, while coefficients B and D model transport effects of, respectively, 
viscosity and heat conduction, and species diffusion. Multi-step reactions may be 
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modeled by the same equations with vectorial reaction variable z G IR m , and coef- 
ficients g, .D, <p, k modified accordingly; see Section 9. For further discussion, see, 
e.g., [CF,FD,GS.l-2,Z.3,Ly]. 

Under different conditions at x = ±00, there can result a variety of types of waves 
solving (1.1): nonreactive gas-dynamical shock and rarefaction solutions (z = 0, or 
z = constant and <p = 0), and traveling combustion waves consisting of weak 
and strong detonations, weak and strong deflagrations, and Chapman- Jouget det- 
onation and deflagration waves, which are limiting cases dividing weak and strong 
branches. Roughly speaking, detonations are compressive waves analogous to shock 
waves in nonreactive gas dynamics, while deflagrations are expansive solutions anal- 
ogous to rarefactions; for a fixed left-hand state, there are weak and strong branches 
of right-hand sides corresponding to waves of each type. Chapman- Jouget waves 
occur at the special parameters for which strong and weak branches coalesce. We re- 
fer the reader to [CF,FD,G,M.4,GS.l-2] for a detailed discussion of these solutions 
of the traveling wave ODE for (1.1) and their roles in Riemann solutions/time- 
asymptotic behavior for the initial value problem under various assumptions on B, 
D. 

Similarly as in the case of "real" , e.g., van der Waals, gas dynamics [BE,MeP], the 
multitude of possible such elementary waves leads to a multitude of possible time- 
asymptotic states, and these must be classified according to stability. However, 
the assessment of stability is a complicated undertaking. Up to now, essentially all 
analyses have been carried out for one of three simplified models: (i) the Zeldovich- 
von Neumann-Doering (ZND) model, for which B and D are set identically zero 
in (1.1) [Er.l-6,CF,LS], (ii) the still further simplified Chapman- Jouget (CJ) or 
"square-wave" model, for which B and D are set identically zero and k is taken to 
be infinite, so that combustion waves become surfaces of discontinuity [Er.7,F,etc], 
or (iii) The Majda model, for which u is taken to be a scalar, and B and D are set 
to 1 and (or sometimes 1), respectively [M.4,LLT,LYi,RV], or the Majda-Rosales 
model [CMR,Sz,LY,Li.l-7], a closely-related cousin in which — z x is substituted 
for zt everywhere in (1.1). At one time, it seems to have been believed, based 
on analysis of the ZND case B = D = (see, e.g., [CF]) that weak detonations 
and strong deflagrations were unstable, the other types stable at least in moderate 
parameter ranges. This conjecture on weak detonation is now widely agreed to 
be false when viscosity and other effects are taken into account, see [CMR,Sz,LY], 
or more general discussion in [FD]. However, rigorous analysis of stability for the 
full model (1.1), or comparison with stability for the ZND or CJ approximations, 
remain important open problems. 

The purpose of the present paper is to initiate a larger-scale study of these prob- 
lems by the introduction, in the simple setting of the Majda model, of new Evans 
function techniques developed recently in the study of stability of viscous shock 
profiles (see, e.g., [GZ,ZH,BSZ,ZS,Z.3]). In contrast to the methods of past analy- 
ses of the Majda model, these techniques were designed for the study of systems, 
u E R n , n > 1, so may be applied also in the case of the full, reacting Navier-Stokes 
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equations. As a roadmap for discussions to follow, we point out that the ZND limit 
B = eBq, D = eDo, e — > is equivalent by the spatial rescaling x — > x/e to the 
small-/? limit k — > with B, D held fixed. The CJ, or square-wave limit in the 
viscous setting (B, D fixed) is ambiguous, corresponding to intermediate values of 
k, but large spatial scale; see discussion of [Z.3], Appendix A. 3. The large k limit 
k — > oo (with, necessarily, U{, u 1 — > U- to allow a connection) corresponds roughly 
to the high activation-energy limit described in, e.g., [B], pp. 22-25, which in the 
ZND setting leads to a square-wave approximation theory, and is generally associ- 
ated with instability, and other anomalous behavior; see, e.g., [Er.7,AT,BL,BN,LS]. 

A tool that has proved useful for the study of stability in the related cases 
of van der Waals gas dynamics and multiphase flow [GZ,Z.6] is a one-dimensional 
stability index originally introduced by J. Evans [E.l-4] in the context of nerve axon 
equations, and generalized in various directions in, e.g., [J.1,AGJ,PW,GZ,BSZ]. 
A topological index relating evolutionary (PDE) dynamics of a traveling wave to 
dynamics of the associated traveling wave ODE, the stability index is based on the 
Evans function [E.l-4,AGJ,PW,GZ,BSZ,etc] -D(A), an analytic function playing 
the role of a characteristic function for the linearized operator L about the wave. 
Defined as a Wronskian of solutions of the eigenvalue equation for L decaying at 
plus and minus spatial infinity, the Evans function vanishes at A if and only if there 
exists a solution of the eigenvalue equation decaying at both infinities, i.e., A is an 
eigenvalue. More precisely, zeroes of D agree in both location and multiplicity with 
eigenvalues of L [GJ.1-2]. 

For traveling waves, there is always an eigenvalue at A = 0, corresponding to 
translational invariance of the underlying PDE, which in the simplest setting is 
multiplicity one. In this situation (which will be the case here), D(0) = 0, but 
-D'(O) 7^ 0. Moreover, if the evolution equation is well-posed in the sense that the 
linearized operator about the wave generates a C° semigroup, then (by standard 
resolvent estimates [Pa]) D cannot vanish for A real and sufficiently large. Since 
.D, properly constructed respects complex symmetry, D(X) = D(X), as does the 
eigenvalue equation itself, we have in particular that D is real- valued when restricted 
to the real axis; thus, sgn D(X) has a well-defined limit as A — > +oo along the real 
axis, which we will denote as sgn D(+oo). The stability index is then defined as 

(1.2) T := sgn D'(0)D(+oo). 

Evidently, the stability index detects the parity of the number of real roots on 
the nonnegative real axis, V positive corresponding to even parity, and V negative to 
odd parity. Since complex roots appear in conjugate pairs, this is in fact the parity 
of the number of all roots (real and complex) in the unstable half-plane Re A > 0, 
and thus gives partial information on spectral stability, defined as nonexistence of 
eigenvalues in this region. In particular, T > is seen to be a necessary condition 
for stability. The value of this index comes from the fact that it can be related to 
geometric information about the phase portrait of the traveling-wave ODE, which 
fact ultimately derives from the correspondence at A = between the eigenvalue 
equation and the linearized traveling-wave ODE. 
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To obtain a concrete result, of course, requires information about the existence 
problem, and this can in general be obtained only in simple situations, e.g. scalar 
reaction-diffusion equations, 2x2 parabolic conservation laws, 3x3 conservation 
laws with real viscosity, etc., for which the connection problem is a planar dynam- 
ical system, or else [J.l,AGJ,etc] in some singular limit for which the connection 
problem can be broken into separately computable "fast" and "slow" problems. 

In this paper we show, in the simple context of the Majda model, that the meth- 
ods introduced in [GZ,BSZ] for the study of stability of viscous shock waves, may, 
with slight modifications, be applied also in the study of stability of detonations 
to: (i) construct an analytic Evans function on the set Re A > for the linearized 
operator about the detonation wave, and (ii) in both the strong and weak detona- 
tion cases, compute an expression for V in terms of quantities associated with the 
traveling-wave ODE. In this simple setting, the connection problem is planar, and 
we can in fact do more, obtaining a complete evaluation of V; the result, for both 
weak and strong detonations, is V > 0, consistent with stability. 

This is consistent with prior results of [L,LLT,LYi] for small-amplitude strong 
detonations in the small-g limit, and of [RV] for arbitrary amplitude strong det- 
onations in the small-Zc (ZND) limit. It is also consistent with results on weak 
detonations of the Majda-Rosales model obtained in [Sz] for small-amplitude waves 
with intermediate k and in [LY] for arbitrary amplitude waves in the large-/? limit; 
however, as far as we know, ours is the first analytical result on stability of weak 
detonations for the Majda model. 

These are partial stability results in that they do not rule out instability; on the 
other hand, they are of general applicability, and are obtained using a relatively 
small amount of information about the system under study. In particular, they 
generalize to the full reactive Navier-Stokes case (see discussion below), whereas 
analyses of the integro-differential Majda-Rosales model clearly do not. Moreover, 
they do have the interesting implication that transition to instability, if it occurs, 
must result from a pair of complex conjugate eigenvalues crossing the imaginary 
axis, typically signaling a Poincare-Hopf bifurcation to a time-periodic solution, 
consistent with the experimentally and numerically observed phenomenon of "gal- 
loping" detonations [MT,FW,MT,AlT,AT,F.l-2, p. 161,BMR,S,Li.6]. It would be 
very interesting to search numerically for such instabilities in the high activation- 
energy limit k — > oo; see, e.g., [Br. 1-2, Br Z] for an efficient numerical algorithm. 

In the course of our development, we also discuss existence of deflagrations for 
a modified "bump-type" ignition function, and the appearance in the ZND limit of 
degenerate, spatially subalgebraically decaying families of weak and strong detona- 
tions, as described for the full equations in [GS.l]. Stability of these waves, and their 
significance in Riemann solutions, are discussed in Sections 8 and 4, respectively. 

Discussion and open problems. In [Ly,LyZ], the methods of this paper are 
extended to the full Navier-Stokes equations of reacting flow, with possibly multi- 
species reaction. In this case, the connection problem is high-dimensional, and the 
actual evaluation of the geometric stability conditions is carried out analytically 
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only in the small-/c (ZND) limit, for which the connection problem has been closely 
studied in [GS.1-2]. However, the geometric condition itself is valid for arbitrary 
model parameters, and could be studied numerically at the same time as the con- 
nection problem; thus, the results again reduce the question of stability to state of 
the art existence theory. Extensions to multi-dimensional stability, in the spirit of 
viscous shock calculations of [ZS], are given in [Z.3,JLy.l-2]. 

As pointed out in [Z.3], Appendix A. 3, it is relatively straightforward using the 
methods of [ZH,Z.2-4,MZ.l-3] to show that spectral stability of nondegenerate det- 
onation waves implies linearized and nonlinear orbital stability, both for the Majda 
model 1 and for the full, reacting Navier-Stokes equations; this will be a topic of 
future work. A very interesting open problem is to carry out a complete spectral 
stability analysis in the spirit of the matched asymptotic analysis of the Majda 
model in [RV] for the full reactive Navier-Stokes equations in the ZND limit, both 
in one- and multi-dimensions. The Majda model has the simplifying feature of 
scalar kinetics; recent singular perturbation-type techniques developed in the shock 
wave context in [PZ,FreSz] may be helpful in attacking the full, system case. As 
discussed in Section 8, it would also be interesting to study further the stability of 
(generically appearing) degenerate weak deflagrations and of degenerate weak and 
strong detonations appearing in the ZND limit. 

In actual detonations, multidimensional, geometric effects become important; for 
example, a converging detonation wave is destabilized, while diverging, or expansive 
waves are stabilized, and this is another important feature to understand; see, e.g. 
[MR,B,BM,M,Li.7]. There is also the problem of trying to understand the bifurca- 
tion to quasi-steady behaviors, such as time-oscillatory "galloping", or "spinning" 
detonations in the multidimensional case, in terms of spectral information given by 
the Evans function. The challenge of this problem, similarly as the original stability 
analysis, comes from the fact that there is no spectral gap between neutral point 
spectrum and essential spectrum of the linearized operator about the wave, so that 
tbifurcations are of mathematically nonstandard type; see the appendix of [BMR] 
for an interesting related discussion in the context of the ZND model. Useful survey 
of the latter two topics may be found in [FD,LS,S,Li.6]. Finally, we mention the 
(presumably numerical) problem of cataloguing possible Riemann solutions within 
the class of viscous profiles in the more general situations discussed in Section 9 of 
multi-species reactions or reaction-dependent equation of state; this appears to be 
an interesting and physically important direction for further study. 

Plan of the paper. In Section 2, we give a brief description of the Majda model, 
and the associated Chapman- Jouget diagram. In Sections 3-4, we discuss existence 
of profiles and the implications for solutions of the Riemann problem. In Sections 
5-7 we construct an Evans function for the linearized operator about the wave, and 
carry out the described analysis of the stability index. In Section 8, we discuss the 
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case of degenerate, spatially subalgebraically decaying profiles and in Finally, in 
Section 9, we describe extensions to multi-species reactions and reaction-dependent 
equation of state. 

Section 2. The Majda model 

Hereafter, we restrict attention to the Majda model 



(2.1) 

or, equivalent ly, 
(2.2) 



u t + f(u) x = u xx + kqip(u)z, 
z t = -kip(u)z, 



(u + qz) t + f(u) x = u«, 
z t = -kcp(u)z, 



a simplified, scalar version of (1.1) simulating the dynamics of one-dimensional 
combustion within a single characteristic family [M.4]. Here, u G M 1 is a lumped 
gas-dynamical variable combining aspects of specific volume, particle velocity and 
temperature, and z corresponds to mass fraction of reactant as before. 
Following [M.4], we take 

(2.3) /'(«)> 0, /»>0, 

and q > 0, corresponding to an exothermic reaction. However, in place of the 
"step-type" function of [M.4], we take a modified, "bump-type" ignition function 



(2.4) ^eC\ | 



<p(u) = for u < Ui or u > u l , 
<p(u) > for Ui < u < u\ 



This choice is motivated by the physical parametrization of temperature with re- 
spect to velocity u in the traveling-wave phase portrait of the ZND model, and 
ensures that the traveling-wave equations for (2.1) agree with the reduced system 
obtained in [GS.l] by asymptotic analysis in the ZND limit; see [GS.l], pp. 979-981. 
In particular, it allows for existence of weak deflagration profiles, as the step-type 
ignition function does not; see Section 3, below. 

More precisely, interpreting u as particle velocity we take 

(2.5) V = ifj(T(u)), 

where T{u) denotes temperature, ip is a standard, step-type ignition function, and 
ip = below some ignition temperature Tj and positive above, and T(u) is qua- 
dratic, concave-down, with T(ttj) = T(u l ) = T{. This agrees qualitatively with the 
physical dependence of temperature on velocity along the one-dimensional flow of 
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the traveling wave ODE for the ZND model. For example, for an ideal gas with gas 
constant 7 > 1, the physical dependence is 

(2.6) T(u) = -^M 2 u 2 + ("fM 2 + l)u, 

where M = \u + — s\/\c+\, < M < 1, denotes Mach number of the specific 
detonation under consideration, c + the sound speed at the righthand state; see 
[LyZ], eq. (1.17). 

We seek traveling waves (w, z) = (u(x — st), z(x — st) connecting states (ii_, Z-) 
and (it.)-, z+ ), with 

(2.7) z + := z(+oo) = 1, z- := z(-oo) = 0, 
and 

(2.8) Ui < U- := u(—oo) < u l ; u+ := u(+oo) < Ui or > u l , 
hence 

(2.9) if- := <p(u-) > 0; </?+:= <fi(u+) = : 

that is, combustion waves moving from left to right, leaving completely burned gas 
in their wake. 

The Traveling-wave equation. The traveling wave ODE for (2.2) is 

^ 2 (u' = f(u) - f(u-) - sqz - s(u - u-), 

\ z' = (k/s)(f(u)z. 

Thus, necessary conditions for existence of a profile are the modified Rankine- 
Hugoniot condition 

(2.11) [/]=«([«] + 9), 
along with 

(2.12) <p + : <fi(u+) = 0, i.e. u + < Ui or > u 1 , 

together assuring that (u + ,z + ) = (u + , 1) is a rest state for (2.10). 

Types of Waves. Examining Figure 1, cases A and B, below, we find that the 
possible solutions of (2.11) may be described as follows. 

Proposition 2.1. For fixed u+, s > s*(-u+), there exist two states U- > u+ for 
which (2.11) is satisfied. For s = s*, there exists one solution. For s < s*, there 
exist no solutions U- > u + . Here, s* > f'(u + ). 
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Figure 1a. Rest states it_ > u + 
(detonation case). 



Figure 1b. Rest states it_ < u + 
(deflagration case). 



Proposition 2.2. For fixed u + , s < s*(u + ), there exist two states U- < u + for 
which (2.11) is satisfied. For s = s* , there exists one solution. For s > s* , there 
exist no solutions U- < u+. Here, s* < f'(u+). 

Figure 1, known as a Chapman- Jouget diagram, is the basis for the standard 
classification of combustion waves. Denote 



(2.13) 



Definition: A combustion wave (u,z) is called a detonation if U- > u+. It is 
called a strong detonation if the Lax characteristic condition holds, 



(2.14) 



a_ > s > a + . 



This corresponds to the larger of the two solutions in Proposition 2.1. It is called a 

weak detonation if 



(2.15) 



s > a_, a + , 



corresponding to the smaller of the two solutions. This type is undercompressive. 
The boundary case 



(2.16) 



a_ = s > a + 



is called a Chapman-Jouget detonation, and has special significance in the the- 
ory: specifically, in idealized circumstances, it is the wave expected to be time- 
asymptotically selected in the "ignition problem" of initial data consisting of a 
large initializing pulse; see [FD,Li.2-4], or discussion at the end of this section. 
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Definition: A combustion wave (u,z) is called a deflagration if U- < u + . It is 
called a weak deflagration if 

(2.17) a+,a_>s, 

corresponding to the larger of the two solutions, and strong deflagrations if 

(2.18) a + > s > a_, 

corresponding to the smaller. The former is again undercompressive, the latter of 
"reverse-Lax" type. The boundary case 

(2.19) a+>s = a- 

is called a Chapman-Jouget deflagration. 

The minimum detonation speed s* and the maximum deflagration speed s* are 
called Chapman-Jouget (CJ) speeds. Note that they determine a minimum strength 
for combustion waves connecting to u+. Indeed, there is a band of speeds (s*, s*) 
about the characteristic speed f'(u+) for which neither detonation nor deflagration 
connections exist, in sharp contrast to the inert-gas, shock wave case. Thus, for fixed 
q > 0, all combustion waves are "strong" in the sense that they differ appreciably 
from acoustic signals. 

Section 3. Existence of profiles. 

The Chapman-Jouget classification scheme concerns possible endstates that may 
be connected by a traveling wave. We next consider the question of existence of a 
traveling- wave profile, reviewing and extending results of [M.4,Sz,Li.3,LY]. 

Linearized rest points. Linearizing (2.10) about the critical points (u±,z±), 
we obtain 

(3 1) ( U \ - ( a± ~ sq ^ ( U 

where 

(3.2) a± := a± — s. 

(Note: because of the structure of (2.10), we cannot without loss of generality 
set s = as in the conservative, shock case.) Here, we have used the fact that 
dip(u+) = Z- = 0, so that dtp(u±)z± = 0. Thus, the eigenvalues associated with 
the linearized equations are 



(3.3) 



7± = a±, (k/s)cp(u±). 
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There are two distinct cases: namely, the generic case that u+ differs from ignition 
temperatures Ui and u l , i.e., the inequalities in (2.9), (2.12) are strict, and the 
degenerate case that u+ lies precisely at ignition temperature U{ or u l . 

In the generic case, referring to (2.9), (2.12), we find that (u+,z+) is a saddle- 
attractor in the detonation case, and a saddle-repellor in the deflagration case, with 
center manifold in both cases lying in direction (qs, a+) tangent to the null-cline 
u' = (recall, <p + = 0), while (u-,Z-) is a repellor in the strong detonation or 
weak deflagration saddle in the weak detonation or strong deflagration case. 

In the exceptional, Chapman- Jouget case ct_ = that weak and strong values for 
U- coalesce, (u-,z-) becomes a saddle-repellor, with center manifold lying in the 
gas-dynamical direction (1, 0); see (3.1). Considering the traveling wave equation as 
a scalar ODE vl = f"(u-)u 2 + 0(\u\ 3 ) — sqz, forced by the solution z = zoe^/ 8 ^ of 
the decoupled exponential growth equation z' = (k/s)z, we find that solutions are 
generically governed by the dominant part u' = f"(u-)u 2 , growing algebraically as 
c/\x\; the sole exception is the unique solution lying along the unstable manifold, 
—sqz ~ (k/s)u, for which both u and z exhibit exponential growth with rate e^^^ 1 . 

Finally, note that the center manifold (CM) at (u+,z+) is a manifold of rest 
points, hence there are no orbits connecting to (u + , z + ) along the CM. That is, for 
u < Ui or u > u 1 , ODE (2.10) reduces to 

f u ' = f(u) - f(u-) - s(u - u-) - sqz, 
U' = 0, 

a scalar ODE with parameter z, for which u+ is a nondegenerate rest point by the 
fact that f'(u + ) ^ s, a consequence of Propositions 2.1-2.2: an attractor in the 
detonation case (with a unique incoming orbit), a repellor in the deflagration case 
(with no incoming orbit). 

Collecting the above discussion, we have the preliminary observation: 

Lemma 3.1. In the generic situation u + ^ Ui, u + ^ u 1 , detonation profiles , if they 
exist, are unique up to translation, decaying exponentially to (u+,z+) as x — > +oo. 
Moreover, weak and strong detonation profiles decay exponentially to (u-,z-) as 
x — > — oo, while Chapman- Jouget detonation profiles generically decay algebraically, 
as c/\x\. For fixed u + , s, there may exist a weak or a strong detonation, but not 
both. Neither weak nor strong deflagration profiles exist. 

In the degenerate case that u+ = Ui or u+ = u 1 , the center manifold of rest point 
(u+,z+) consists of equilibria only on the "upper" side z > z+; indeed, it is not 
difficult to see that flow along the CM is attracting in either the detonation (u = Ui) 
or deflagration (u = u l ) case, with rate of approach slower than any algebraic order. 
For, approximating z ~ 1 and u — u + ~ (z — z + )(sq/a + ), we may estimate the order 
of approach by consideration of the scalar ODE 



(z - z+)' = (k/s)(p(u+ + (sq/a + )(z - z+)) > 0. 
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In the detonation case a + < 0, for example, this becomes z' = (k/s)(p(ui + 
(sq/\a+\)\z\), where z := z — z+. Noting that <f(ui + h) = o(h k ) for any alge- 
braic order k, we obtain the result by comparison with ODE z' = \z\ k for z < 0. In 
particular, deflagrations now become possible in principle, while detonations now 
become possibly nonunique. More precisely, we have: 

Lemma 3.2. In the degenerate case u + = ui strong detonation orbits, if they exist, 
occur as a one-parameter family of subalgebraically decaying orbits bounded on the 
upper side by a unique exponentially decaying strong detonation, and on the lower 
side by a pair of orbits consisting of a gas-dynamical shock orbit followed by a sub- 
algebraically decaying weak detonation orbit, decaying exponentially to (u+,z+) as 
x — > +oo. Weak detonation profiles, if they exist, are still unique up to translation, 
but may be either subalgebraically decaying as x — > +oo, in which case they bound 
a one-parameter family of strong detonation orbits as described above, or exponen- 
tially decaying as x — > +oo, in which case there is no strong detonation profile 
connecting to u+. In the degenerate case u+ = u l , there always exists a unique, 
weak deflagration profile, decaying subalgebraically as x — > +oo; strong deflagra- 
tion profiles, however, do not exist. The Chapman- J ouget cases are similar, but 
generically exhibit algebraic decay ~ c/\x\ as x — > — oo, as described in Lemma 3.1. 

We defer the proof to the following subsection. Exponential decay is what is 
needed to apply the general Evans function machinery of [GZ]. The "sonic" case for 
shock waves, analogous to the Chapman-Jouget case has been treated in [H. 1-3, HZ] 
by explicit calculation, and it appears likely that this approach should generalize 
to combustion; however, we will not treat that issue here, ignoring for the moment 
the case of Chapman-Jouget detonations in our Evans function analyses. In the 
case that u+ = U{ or u l , for which profiles in general approach (it + , z+) as x — > +oo 
at subalgebraic rate, it is not possible by current techniques to define the stability 
index, and a different approach must be taken in analyzing stability. 

Phase Plane Analysis. A more detailed description of existence properties 
may be obtained by examination of the phase plane, as depicted in Figure 2, below. 
Following [LY], consider the nullcline 

(3.4) F(u) := f(u) - f(u+) - s(u - u+) - sq(z - z+) = 

on which u' = 0. By convexity of /, this intersects the z = line at the two points 
U- and the z = 1 line at the two points u + described in Propositions 2.1 and 2.2; 
from here forward, denote these by u+ < U- < u~ < u + . By positivity of sq, the 
nullcline F = is a graph over u, whence on the physical, invariant region z > 0, 
the region F < above the nullcline is attracting (recall, z' > for z > 0) so 
that any flow entering this region for u > Ui remains there and (since u' = F < 0) 
eventually crosses the vertical half-line u = Ui, z > 0. The value z = z at which 
it strikes at u = ui will be the terminal z- value, since z' = for u < Ui, and the 
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Figure 2. Phase portrait. 

terminal w-value u will be the unique point u < ui at which the nullcline intersects 
z = z. 

^From these basic considerations, we may conclude already that strong deflagra- 
tions are under no circumstances possible, since these require a passage from U- 
to u + , and these points are separated on z > by the region F < from which 
orbits can escape only to the left. Weak deflagrations u~ — > -u + , if they exist, must 
lie entirely within the complement F > 0, hence are monotone increasing in both z 
and u. It is easily seen that they in fact exist if and only if u + = u l (recall, we have 
already shown that they do not exist for u + > u l ). More generally, there always 
exists a connection from (u~,0) to the rest point (u l ,z*) determined by the inter- 
section of the nullcline F = and the vertical line u = u l ' . For, tracing backward 
along the center manifold leading to (u l ,z*), we see that it is trapped in F > 0, 
z > 0, hence must terminate at the corner (u~, 0) of this wedge. 

Likewise, weak detonations U- — > u + , if they exist, are monotone in both z 
(increasing) and u (decreasing), lying entirely within F < 0. For, the unique orbit 
originating from saddle U- into z > lies along the unstable manifold, pointing 
into the region F < 0, from which it can exit only at a rest point along the left hand 
edge of the nullcline F = 0. Moreover, denoting by z the value of z at which this 
special orbit strikes Ui, we find that: (i) there exists a weak detonation profile if 
and only if z = 1, necessarily unique, (ii) there exists a unique strong detonation 
profile if and only if z < 1. (iii) there exists a one-parameter family of degenerate 
strong detonation profiles if and only if u+ = Ui, z = 1, and the weak detonation 
profile approaches (u+, 1) along the center manifold; if the weak detonation profile 
approaches along the stable manifold, then there are no other connections. (Note 
that z > 1 when u+ = Ui, since the nullcline is a lower barrier for the unstable 
manifold.) (iv) for z > 1, there exist no detonation profiles, neither weak nor 
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strong. 

For, in case (iv), this orbit (the unstable manifold from (tt_,0)) serves as a 
barrier separating (u+,1) from all orbits in z > originating from (tt + ,0). This 
is also the case when there exists a weak detonation profile connecting to (u+, 1) 
along the stable manifold. If z < 1 and the unstable manifold of (w_,0) does not 
coincide with the stable manifold of ("u_|_, 1), then the latter must lie entirely above 
the former; tracing it backwards, we find that it must either approach the rest point 
(u~ , 0) within F < 0, corresponding with a monotone strong detonation profile, or 
else exit F < along its righthand edge, in which case they must approach (w~, 0) 
from within the wedge F > 0, z > 0, by the same argument used to prove existence 
of weak deflagrations. If, along with this nondegenerate strong detonation profile, 
there exists also a degenerate weak detonation profile entering (u+, 1) along a center 
manifold, i.e., u+ = Ui and 5 = 1, then trapped between these two orbits is a one- 
parameter family of degenerate strong detonation profiles also entering along the 
center manifold. This verifies claims (ii)-(iv), while claim (i) is evident. 

We have now verified the claims of Lemma 3.2; more, we have characterized 
the existence of weak and strong detonation profiles in terms of the height z at 
which the unstable manifold of (w_,0) strikes Ui (and in degenerate cases also the 
direction): equivalently, in terms of the Melnikov separation function 

(3.5) d(u + , s, k, q) := z — 1. 

The following monotonicity properties will be helpful in completing our description 
of existence. 

Proposition 3.3.. The function d(-), more generally, the height of the unstable 
manifold of (w_,0) at any value Ui < u < U-, is monotone nondecreasing in u + , 
k, and, for z < 1, in q, and is nonincreasing in s. Moreoever, this monotonicity is 
strict except in the special case that (ui,z) is a rest point, i.e., lies on the nullcline 
F = 0. 

Proof. To establish monotonicity, we show in each case that the vector field 
(F, (7), G := (k/s)<f(u)z, determining the flow of the traveling wave ODE, within 
the region F < 0, z > 0, Ui < u < U- rotates clockwise strictly monotonically with 
respect to the parameter j3 under consideration, or, equivalently, 

(3.6) **==**(g 9GmO S °' 

with < corresponding to monotone increase, and > to monotone decrease. For 
example, 5k is simply FG/k < 0, 5 q is Gs(z — z + ) < 0, and 5 U+ = Ga + < 0, by the 
fact that a+ < (weak detonation). Finally, 

l -det^ " 5(Z 7w + L"SV" +) ) (V*)(F -(/(«)-/(«+)) \ 

\G -(k/s 2 )ip(u)z J \G J 

= (-G/s)(2F -(/(«)-/(«+)) 

> 0. 
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This gives strict comparison principles on ui < u < U-, not only for the unstable 
manifolds, but for any orbits lying in F < 0, z > (1 > z > in the case of q), 
Ui < u < u l . Noting that du-/ds < 0, du-/dq, du-/du + > 0, and du-/dk = 0, 
we thus obtain the claimed strict monotonicity on Ui < u < U-. It follows that 
there is an orbit lying strictly between two stable manifolds under consideration for 
Ui < u < and this implies strict monotonicity at u = Ui unless this third orbit 
collides with the stable manifold at Ui, which can only happen if (m,z*) is a rest 
point. | 

(Note: the above argument applies also in the Chap man- Jouget case.) 

Remark 3.4. In the nondegenerate case u + ^ Ui, we may express the derivatives 
of d at a profile d = in the standard way (see, e.g., [GH,HK]) as Melnikov integrals 
(3.8) 

dd/OP = -J°° e-fo ^ (8(F,G)/8(u, z )K*, z )(y))dy det ^ fj ^ (u,z)(x)dx, 

to find by the same calculation that the derivatives dd/ d/3 have strict signs as well: 
dd/du +: dd/dk, and dd/dq > 0, while dd/ds < 0. In particular, we have 



(3.9) 



dd/ds= [ X e-SZfr+WWvMdetf** & ~ u +] + «> " ) (x)dx 

/oo 
e -/ *(a+M(^))^(^/ s )(2u' - (f(u) - f(u+)))(x) dx < 0. 
-oo 



Remark 3.5. Alternatively, following [M.4], we may fix the lefthand state state 
U- and consider d as a function of (k,q,U-,s). It is a straightforward exercise to 
verify that d then becomes monotone nondecreasing in U- and k, but nonincreasing 
in q (now for all values of z). Monotonicity in s is lost. 

Our main conclusions regarding existence are summarized in the following two 
propositions. 

Proposition 3.6. For each u+, there holds one of two possibilities: either (i) 
for each speed s greater than or equal to the minimal, Chapman- Jouget detonation 
speed s* for which detonations can occur, there exists a single exponentially decaying 
strong detonation profile, or (ii) there exists a threshold speed s > s* at which there 
exists a (necessarily) unique exponentially decaying weak detonation profile, above 
which there exists a single exponentially decaying strong detonation profile, and 
below which there exists neither weak nor strong profile. Where it is defined, s* is 
monotone increasing in u + . If u + ^ Ui, then strong detonation profiles, when they 
exist, are the only detonation profiles that occur. If u + = Ui, on the other hand, 
then whenever there exists an exponentially decaying strong detonation, it bounds 
from above a degenerate family of strong detonation profiles, bounded below by a 
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degenerate weak detonation profile and a gas-dynamical shock profile, as described 
in Lemma 3.2. Weak detonations are always monotone in u (decreasing) and z 
(increasing) . 

Weak deflagrations exist if and only if u + = u 1 , and are monotone increasing in 
both u and z; moreover, they are always degenerate, decaying sub-algebraically as 
x — > +00. Strong deflagrations do not occur. 

Proposition 3.7. For fixed model parameters f , <p, q, there exists ko > such 
that, in Proposition 3.6, only case (i) occurs for < k < ko, while, for k > ko, case 
(ii) always occurs for some choice of (u+, s). 

Proof of Proposition 3.6. It is straightforward to see, in the infinite-speed 
limit U- — > u + + q, that z < 1 if u + 7^ Ui, and z = 1 if u+ = Ui, with the weak 
detonation profile entering u+ = ui along the center manifold. If u + + q < Ui is out 
of range, then we consider instead the maximum speed limit U- — > Ui, observing 
that z — > is again less than 1. As speed s is decreased, z increases monotonically, 
so that either z < 1 for all s, with the unstable manifold of U- lying always below 
the stable manifold of u+, or else there is a unique transition point s* as described 
in the theorem at which the two manifolds coincide, below which z > 1 and there 
are no detonation connections. That is monotone in u + follows from the fact 
that z decreases with s but increases with u + . The remaining conclusions about 
detonations, and the conclusions about weak and strong deflagrations, follow from 
the discussion above Proposition 3.3. B 

Proof of Proposition 3.7. As k — > 0, the unstable manifold of (w_,0) ap- 
proaches the nullcline F = 0, as can be seen either by direct calculation, or using 
singular perturbation techniques as in [GS.l]. Likewise, the stable manifold of rest 
point (v,i,z*) (recall: the intersection of the nullcline and u = Ui) approaches a 
horizontal line. Thus, the unstable manifold of U- is trapped between the nullcline 
and the stable manifold of (ti;,2:*), hence must approach (ui, z*) along the center 
manifold tangent to the nullcline, decaying subalgebraically as x — > +00; see the 
discussion of case u+ = Ui for details. (In particular, it connects to (u+, 1) if and 
only if u+ = u^.) ^From these observations, and z = z* < 1, we find that we are in 
case (i). Moreover, this argument is uniform in model parameters, yielding a global 
result for k less than some threshold k as described. That case (ii) occurs for all 
k > ko follows by monotonicity of z with respect to k. | 

Section 4. Riemann solutions and the CJ shift. 

We conclude our discussion of existence by cataloging solutions of Riemann prob- 
lems involving waves with viscous profiles. In a Riemann problem, we prescribe 
data Ul = (ul,zl) for x < and Ur = (ur,zji) for x > 0, and seek a sequence 
of waves (Ul,U\), (Ui, U2), • • • , (U m , Ur) with increasing speeds, progressing from 
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left state Ul to right state Ur, consisting of gas-dynamical shock or rarefaction 
waves, weak or strong detonations, or weak deflagrations, each possessing a viscous 
profile, for which the endstates Uj = (uj, Zj), j = 1, . . . , m, R are individually valid 
time-asymptotic states, i.e., Zj = if itj < Uj < u l . As discussed, e.g., in [AMPZ], 
such solutions represent possible time-asymptotic states for solutions of (2.1), (2.2), 
as obtained by formal, matched asymptotic expansion with the usual, hyperbolic 
scaling x, t — > 

By the one-sided nature of the Majda model, it is easily seen that the only 
interesting situation is the one zl = 0, zr = 1 described in the introduction. For, 
the reverse situation zl = 1, zr = clearly admits no solution, there being no waves 
along which z increases, nor, for the same reason, does the situation zl = zr = 1 
with ul and ur lying in different components of (— oo,itj] U [u l ,+oo). In these 
cases, the asymptotic behavior has a different, diffusive scaling not captured by the 
Riemann solution. For zl = zr = on the other hand, or zl = zr = 1 and ul 
and ur lying in the same component of (— oo, Ui] U [u\ +oo), the solution is just the 
gas-dynamical solution for the conservation law ut + f(u) x = u xx with the value of 
z held fixed. We therefore fix zl = 1, zr = for the remainder of our discussion. 

Case I. (deflagration: ul < u l < u r ). Denote by ucj the Chapman- Jouget 
deflagration speed associated with u + = ur. Then, there are two subcases. Ia. For 
ucj < ul < u l , the solution consists of a weak deflagration from ul to u l , followed 
by a (possibly zero strength) fluid-dynamical rarefaction with z = 1 from u 1 to ur. 
Ib. For ul < ucj, the solution consists of a (possibly zero strength) fluid-dynamical 
rarefaction from ul to ucj, followed by a Chapman-Jouget deflagration from ucj 
to u l , followed by a fluid-dynamical rarefaction with z = 1 from u l to ur. Note 
that weak deflagrations may be followed by gas-dynamical waves, whereas the CJ 
deflagration may be both followed by and preceded by another wave, so long as the 
preceding wave is a rarefaction. 

Case II. (standard detonation: ul < u l , u r < Ui). In this case, the solution struc- 
ture depends on whether we are in case (i) or case (ii) as described in Proposition 3.6. 
Denote by ucj the Chapman-Jouget detonation speed associated with u + = ur. 
If we are in case (i), then there are two subcases. Ha. For ucj < ul < u 1 , the 
solution consists simply of a strong detonation from ul to ur. lib. For ul < ucj, 
the solution consists of a (possibly zero strength) gas-dynamical rarefaction from 
ul to ucj, followed by a Chapman-Jouget detonation from ucj to ur. 

If we are in case (ii), there are again two subcases, but with the special role of 
the CJ detonation, which no longer has a profile, now played by the unique weak 
detonation possessing a (nondegenerate) profile. Denote by u* < ucj the special 
left state U- that is connected to ur by a nondegenerate weak detonation profile 
with speed s*, and u* > ucj the corresponding value of u~: i.e., the unique state 
that is connected to by a gas-dynamical shock with the same speed s*. IIa\ For 
u* < ul < u 1 , the solution consists simply of a strong detonation from ul to ur. 
lib'. For ul < u* , the solution consists of a (possibly zero strength) gas-dynamical 
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shock or rarefaction, according as ul is > or < -U*, from ul to u*, followed by a 
weak detonation from u* to ur. 

Case III. (limiting detonation: ul < u % , u r = u 1 ). In this degenerate case, there 
exist, along with each nondegenerate strong detonation profile from u~ to Ui, a 
degenerate weak detonation profile from U- to Ui. This leads to nonuniqueness of 
Riemann solutions for ul on the interval between uqj or u* (in case (i) or (ii), 
respectively) and u l . Namely, along with the solutions described in case II, we 
also have solutions obtained by substituting for any strong detonation with speed 
s > s* from u~ to u+ a gas-dynamical shock with speed < s from u~ to u < w_, 
followed by a degenerate weak detonation with speed s > s > s* from u to u + 
(recall: both strong detonation and degenerate weak detonation profiles exist for 
any speed > s*). There is a milder nonuniqueness at the level of profiles, since any 
nondegenerate strong detonation profile may be replaced by a degenerate one with 
the same endstates. 

Case IV. (ul > u 1 , ur > Ui). It is easily deduced that this problem has no 
solution. For, in order that z increase from to 1, there must be an intermediate 
state within (ttj, u l ), and this can only be reached from ul by a single gas-dynamical 
shock. The only wave that can follow a shock is a weak detonation, and no wave 
can follow a detonation. But, ur, as the right endstate of a weak detonation must 
then satisfy ur < Ui, a contradiction. 

Case V. (ul > u l , ur < Ui). By the discussion of the previous case, the only 
possible solution is a gas-dynamical shock from ul to u±, followed by a weak det- 
onation from u- = ui to u+ = ur having the special property u + < Ui < it_ < 
u 1 < ul < u~ . This never occurs in detonation case (i), but may occur in case (ii) 
within certain parameter range. 

Note that the solvable cases I III contain the basic propagation problem described 
in the introduction, of a combustion wave moving from a burned, superignition 
state at — oo into an unburned, subignition state at +oo. Likewise, cases II-III 
and IV contain the ignition problem Ul = (ul,0), Ur = (ur,1), where ul,ur < 
Ui or ul.ur > u l : a one-sided version obtained by left-right symmetry of the 
ignition problem Ul = 1) = Ur for the full, reacting Navier-Stokes equations, 
corresponding to a large, pulse-type excitation, initiating combustion, of a quiescent 
background state. We do not have a useful interpretation of case V. In cases I II, 
the Riemann solution is uniquely determined within the class of viscous profiles, 
while in case III it is not. If we restrict to the class of nondegenerate (exponentially 
decaying) detonation profiles, then the Riemann solution is uniquely determined 
also in case III. A special role in the solution structure is played by the detonation 
profile with minimum speed: the CJ detonation, in case (i), the weak detonation 
profile with speed s#, in case (ii). 

Unique solvability in the class of exponentially decaying profiles suggests that 
nondegenerate detonation profiles, both weak and strong, at least in usual circum- 
stances are stable. Change in stability would presumably signal bifurcation to more 
complicated time-asymptotic behavior, perhaps involving time-oscillatory "gallop- 
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ing" solutions. 

Remark 4.1. (CJ shift) In both propagation and ignition Riemann problems 
considered above, appearance of a weak detonation profile, case (ii), is associated 
with a shift in the solution structure from CJ to weak detonation, in contrast with 
the behaviour predicted by the ZND model. Such a shift is indeed observed in 
experiments, for appropriate parameter regimes [FD]. Proposition 3.7 states that 
no such shift occurs for k < k : that is, it validates the ZND picture for k merely 
small, and not only in the k — > limit. This important observation was made first 
in [GS.l], in the larger context of the full reacting Navier-Stokes equations, and 
answers in the negative a conjecture of Majda [M.4]. 2 

Remark 4.2. As should be apparent from the above discussion, it is monotonic- 
ity of d with respect to s that leads to unique solvability of the Riemann problem 
in the presence of weak detonation profiles (detonation case (ii)). Likewise, a local 
analysis in the spirit of [ScSh,ZPM,Fre.l,GZ,ZS] shows that dd/ds ^ is equivalent 
to linearized well-posedness of the Riemann problem about weak detonation data. 



Section 5. Construction of the Evans function. 

We now turn to the question of stability, beginning by a careful construction of 
the Evans function. Consider a nondegenerate, i.e., spatially exponentially decaying 
traveling wave profile (u(x — st) of any type, satisfying the nonsonicity assumption 
a± 7^ s. The linearized equations of (2.1), (2.2) about u(x — st), in moving coordi- 
nates x = x — st, are 

^ | u t - q{kip'{u)uz - kip{u)z) + {au) x = u xx 

\ z t — sz x = — k(p'(u)uz — k(p(u)z, 



(5.2) 



u t + qz t + (au) x - sqz x = u a 
zt — sz x = — k(p'{u)uz — k(p(u)z, 



where 

(5.3) a := a(u) — s, a(u) = df(u). 

The associated eigenvalue equation is thus 



(5.4) 



u" = (an)' + Xu — q(k(p' (u)uz + k(p(u)z), 
z' = (l/s)(k(p'(u)uz + k<f(u)z + Xz), 



2 In the Remark below Theorem 1 of the reference, the conjecture that q^ R > q for any kg > 0: 
in our notation, that z > 1 for u+ = Uj, for any k > 0, which would imply the existence of 
nondegenerate weak detonation profiles for u+ less than but sufficiently close to Uj. 
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(5.5) 



(5.6) : 



or, alternatively, 

u" = (an)' + Xu + q(Xz — sz'), 

z' = (l/s)(k<f'(u)uz + k(p{u)z + Xz), 

and the limiting systems at x = ±oo are 

u" = a±u' + Xu — qktp±z, 

z ' = (l/ s )(k(f±z + Xz). 

(Recall: <p- > 0, Z- = dtp+ = <p+ = 0.) Here, a± ^ 0, by the nonsonicity 
assumption a± ^ s. 

Seeking solutions u = e^U, z = e^ x Z of (5.6)-|-, we obtain the characteristic 
equation 

> 2 -^«±-A -qk<p± \ fU\ _ /0 N 

a* - (k/s)(p± - X/s J \ZJ ~ [p, 

This is readily solved using its block triangular form, to yield growth rates (eigen- 
values) 



(5.7) 



± 



a± ± Ja 2 ±+ AX 

(5.8) ± fi= V — , (k/s)<p±+\/s, 

>o 

associated, for all except finitely many points A where eigenvalues may coincide, 
with respective normal modes of form 

fu\ ( U \ ( 1 \ /* 

(5.9) ± \u'\ = e^ x I nU 1 = e^ ±x I /j± 

Inspection of formulae (5.8)±-(5.9)± yields 

Lemma 5.1. On {Re A > 0}, there are two unstable and one stable modes of 
(5.7)± at either of ±oo, 

fif < < (i£ t fj,f; 

moreover, there exist analytic choices of solutions of form 

u\ / u \ / * 

u> \ (A) = 

spanning the associated stable/unstable manifolds. Moreover, each of these analytic 
functions has the complex symmetry f(z) = f(z). 

Proof. The first claim is immediate from the formulae; the second then follows 
as described in [GZ,Z.3] by a standard lemma of Kato [K], which asserts that ana- 
lytic subspaces on a simply connected domain possess an analytic choice of bases. 
The third statement follows exactly as in [GZ,Z.3] from the observation that the 
construction of Kato preserves complex symmetry. a 




Further, we have 
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Lemma 5.2. The functions n- , u- and z- can be analytically extended onto a 
region 

(5.10) tt:= {Re \> -e 1 -0 2 \lm\\ 2 }, 9 U 6 2 >0. 

Moreover, are analytic on a neighborhood of X = 0, with associated (distinct) 
analytic eigenvectors. 

Proof. A closer look reveals that Lemma 5.1 holds true on O \ 5(0, r) for 
arbitrary r > and 6j(r) > sufficiently small. Thus, it is sufficient to show that 
there exists an analytic extension on B(0,r): in particular, to verify the second 
claim. At x = — oo, > 0, and so, at A = 0, the reactive root is ji = (k/s) +X/s > 
0, while the kinematic roots are and a + ^ 0; thus, there is a spectral gap between 
the single zero mode and strong stable or unstable modes, and so each of these 
may be analytically continued through the origin. At x = +oo on the other hand, 
<p + = 0, and (5.7) ± becomes block diagonal, decoupling into kinematic and reactive 
blocks. Within blocks, the roots ji are distinct at A = 0, and the result follows as 
before. | 

Remark. For the purposes of this paper, it is sufficient to establish Lemma 5.2 
on the smaller set O = {A : Re A > 0}. 

Applying the Gap lemma of [GZ,KS], 3 we obtain, finally, 

Corollary 5.3. There exist solutions (A, x) , z^(X,x) of eigenvalue equation (5.5), 
analytic on O, (5.10). Moreover, on {Re A > 0}, (u^,z^) spans the stable manifold 
at +oo, and {(u 2 ,z 2 ), ( u 3~, z 3~)} the unstable manifold at — oo. 

Definition: Following the development of [GZ] we define the Evans function as 

U+' «2 ' U s' 

Z l Z 2 H I \ x=0 

Proposition 5.4. The Evans function D is analytic on O; on the subdomain 
Re A > 0, its zeroes correspond precisely to eigenvalues of the linearized opera- 
tor L. Moreover, it possesses complex symmetry D(z) = D(z); in particular, D(X) 
is real for real X. 

Proof. Analyticity and complex symmetry are inherited from the properties of 
the columns of the matrix from which D is evaluated. Recalling that, on Re A > 0, 
the columns of that matrix span the decaying manifolds of the eigenvalue equation at 
±oo, we find, further, that vanishing of D is equivalent to nontrivial intersection of 
these manifolds, i.e., existence of a solution decaying at both spatial infinities. Thus, 
zeroes of D correspond in location with eigenvalues of L. That they correspond also 
in multiplicity follows from a more detailed calculation of Gardner and Jones [GJ.l- 
2]; see also Section 6 of [ZH]. B 



Here, we make use of the assumed spatial decay of the wave, a hypothesis of the Gap lemma. 



GREGORY LYNG AND KEVIN ZUMBRUN 



21 



Corollary 5.5. A necessary condition for linearized stability is that the number of 
zeroes of D in the unstable half-plane {A : Re A > 0} be zero; in particular, that it 
be even. 



Section 6. Stability analysis for strong detonations. 

Specializing now to the strong detonation case, ct_ > > aq_, we complete our 
stability analysis by a computation of the stability index V := sgn D' (0)£)(+oo). 

Analysis at A = 0. At A = 0, (5.5) reduces to the linearized traveling wave 
ODE 

, u" = (au)' — sqz' 

(6.1) 



z' = (l/s)(k(p'zu + k(f{u)z). 
Take without loss of generality 

(6.2) (u^,z^) = (u^,z^) = (u x ,z x ), 

both exponentially decaying modes. The solution (u^ , z 2 ) is likewise asymptotically 
decaying at — oo, 

(6.3) u 2 (~ °°) = z 2 (~ °°) = 0, 

since, for a_ > 0, the unstable modes at — oo remain strictly unstable as A — > 0. 

Integrating the first equation of (6.1) from ±oo to x, we thus obtain the homo- 
geneous first-order equation 

(6.4) ± u' — au + sqz = 0, 

satisfied by each of (-u^,^), (itj,2^"), (u^,z^). 
Proposition 6.1. D(0) = 0, while 

D'(0) = 7A 

where 

(6.5) 7= lim det( U2 _ Ux V-M)e^-^ (a(a:)+(fe/s) ^ (a:)) . 

M^+oo \z 2 z x J 

and 



(6.6) A = [it] + g. 
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Proof. -D(O) = is immediate from (6.2). Using the Leibnitz rule, we obtain as 
usual 



(6.7) D'(0) = det 

where y := y~ — y + , a := a~ — a + , a 
(d/d\)(u^, z$), satisfy the variational equations (obtained by differentiating (5.5)) 

(6.8) 



u 2 


V 






u 2 


y' 


Z 2 


a 


<T + ) 





y ±" = (ay ± ) / - sqa ±f + u x + qz x , 
a ±; = (l/s)(kip'(u)y ± z + k(p(u)a ± + z x ), 



with 

(6.9) (y ± ,a ± )(±oo) = 0. 
Thus 

(6.10) y ±! — ay ± + sqa ± = u — u± + q(z — z±), 
and 

(6.11) y - ay + sqa = [u] + q, 

where [u] := (u + — it_). Using (6.4)±, (6.10) to partially eliminate the 2nd row in 
(6.7), we obtain 

' u x u 2 y 
D'(0) = det I [u] + q 



(6.12) 

= 7 A, 



z x z 2 <J 



-det( "1 ) ([u]+q) 



( U 2 


U x \ | 


\ Z 2 


z x J 



where 7 := det ( "i. _ x ] | x=0 is a Wronskian for the linearization about (u, z) 



of the traveling wave equation (2.10). The evaluation (6.5) then follows by Abel's 
formula.! 



Analysis as A — > +00. Similarly as in [GZ], we next establish 
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Proposition 6.2. As A — > oo afong t/ie rea/ axzs, 

(6.13) sgn D(A) -> sgn w+(+oo) det ( U l U l ) (-oo) 

Proof. Following [GZ], consider the frozen-coefficient version 

( / = a(u)u' + a(u) x u + Xu — qkipz, 

(6.14) < 

t z' = (l/s)(k(p'(u)uz + k(f(u)z + Xz), 

of (5.5) obtained by evaluating u at a fixed point xq, and as usual seek solutions 
u = e^U, z = e F Z: that is, solutions of characteristic equation 

, , f fi 2 - pa- a x - X -qk^p \ fU\ _ (0 

[ } V -(k/s)<p'z fi - (k/s)<p -X/s) \z)~ VO 

It is a straightforward exercise to show that, as A — > +oo along the real axis, 

(6.15) has two kinematic solutions 



fj, 



iA 1 / 2 , 



Z I 



or (u, u', z) ~ e Atx (l, iA 1 / 2 , 0) in phase plane variables, and one reactive solution 

or (u,u',z) ~ e Ma; (0,0, 1) in phase variables. Thus, the stable and unstable sub- 
spaces of the frozen-coefficient equations have a uniform spectral gap of order A 1 / 2 , 
and, carefully applying the tracking lemma of [GZ,ZH,Z.3], we may conclude that 
the stable and unstable manifolds of the variable-coefficient equations (5.5) track 
the corresponding subspaces of the frozen-coefficient equations to angle A -1 / 2 . See 
[Z.3] and especially [MZ.3] for similar calculations. 
^From these considerations, we find that 





(ut 


1 


uf 




u + 






u 2 


u 3 


Z 2 


z 3 




1 

-VX \ T)_{x), 

1 
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where rj + (x) is a real nonvanishing scalar function, and rj-(x) is a real nonsingular 
2x2 matrix function. Thus, 



from which we may deduce the claim by reality of D and nonvanishing of 77+ , det r\- . 

Alternatively, following [BSZ], we may first deduce by by standard Garding-type 
energy estimates that there exist no zeroes of D/eigenvalues of L for sufficiently 
large real A, uniformly, independent of (bounded) model parameters. Performing 
the homotopy /, k — > reduces the eigenvalue equations to the simple case 



for which the claim follows by an elementary direct calculation. a 

The stability index. Combining our calculations at and +oo, we have 

Corollary 6.3. For nondegenerate strong detonations, the stability index V := 
sgn D'(0)D(+oo) is well-defined, and satisfies 



Thus, the number of zeroes of D on the unstable half-plane is even, consistent with 
stability. 

Proof. D(X) is real for real A, and, by the previous lemma, nonvanishing for 
real A sufficiently large. Thus, sgn D(X) is independent of A for real A sufficiently 
large, and sgn D(+oo) is well-defined. Moreover, since D is analytic, the number 
of positive real roots of D is even or odd according as V is positive or negative 
(indeterminate parity if T = 0, corresponding to D'(0) = 0). In fact, since nonreal 
roots occur in conjugate pairs, by complex symmetry of D, the sign of V determines 
the parity of the number of zeroes of D on the entire unstable half-plane {A : Re A > 
0}. Thus, it is sufficient to verify (6.16) in order to establish the claim. 

To establish the formula V = sgn ([u] + q)j 2 u x (+oo), it is sufficient to show 
that expression (6.13) is independent of A for real A > 0; evaluating at A = and 
combining with the formula of Proposition 6.1 for -D'(O), we then obtain the result. 
Here, we are using also the fact that 




= 77+ det r\- 




u" = Xu, 
z' = (X/s)z, 



(6.16) 



r = sgn ([u] + q)-f 2 u x (+oo) > 0. 



sgn 7 = sgn det ( "L i ) (0) = sgn det 

V^2 * 3 / 




( 




) 



(-oo)^O 



GREGORY LYNG AND KEVIN ZUMBRUN 



25 



as a nontrivial (since its columns are by construction independent as x — > — oo) 
Wronskian for the linearized traveling-wave ODE. 

To establish independence of (6.13) with respect to A, we may prove a projec- 
tion lemma as in [GZ], namely that projection onto u coordinate is full rank on 
the unstable subspace, and projection onto (u, z) coordinates is full rank on the 
stable subspace of the coefficient matrix of the limiting eigenvalue equation. This is 
straightforward in the present case, using the limiting structure (1, fx, 0)* of unstable 
eigenvectors, and 



of stable eigenvectors at x — > ±oo. 

Finally, the conclusion sgn ([u] + q)^f 2 u x (+oo) > follows from the facts that: 

(i) u x (+oo) < for all strong detonation profiles, by our earlier phase plane analysis, 

(ii) [u] + q = (l/s)[/(it)] < by [u] < 0, s > 0, and the fact that / is strictly 
monotone increasing, and (iii) 7 ^ 0. | 

Remark 6.4. A similar result has been shown in the ZND limit (k sufficiently 
small) in [Ly,LyZ], for the full, Navier-Stokes equations of reacting flow with an 
ideal gas equation of state. Multi-dimensional analogues are given in [Z.3,JLy.l-2]. 

Remark 6.5. A determinant analogous to the Evans function was used by 
Erpenbeck [Er.2,Er.4] to study stability of discontinuous strong detonation fronts 
within the framework of the ZND model; moreover, he obtained an interesting sta- 
bility criterion by examining the high frequency limit, opposite to what we do (and 
outside the low-frequency regime in which the ZND model is expected to describe 
model (1.1)) [Er.3]; for further discussion, see [Z.3], Appendix A. 3, or [JLy.2]. The 
stability criteria of Erpenbeck-Majda for shock waves [Er.l,M.l-3] , of Erpenbeck for 
detonations [Er.l], and the Evans function criterion of Evans-Alexander-Gardner- 
Jones for reaction-diffusion equations [E. 1-4, J. 1, AG J] are thus seen to be facets of 
the same general principle/stability criterion, applied to equations of varying type. 
For further discussion, see [ZS]. 



We next consider the interesting weak detonation case ct_, a+ < 0, for which 
the speed of the detonation everywhere exceeds the rate of propagation of gas dy- 
namical signals. These correspond to undercompressive, saddle-saddle connections 
and connect u + only to special states it_, as discussed in Section 3. As usual, we 
assume nondegeneracy, i.e., spatially exponential decay. 

Repeating the stability analysis of the previous section, we now find that, at A = 
0, the mode (u^jZ^) is not asymptotically decaying, but asymptotically constant, 




Section 7. Stability analysis for weak detonations. 



(7.1) 
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as x — > — oo (see (5.7)). Thus for j = 2, we obtain in place of (6.4) the equation 
(recall, ct_ < 0) 

(7.2) U2' — 0-U2 + sqz^ = —ct- = \a-\. 
Proposition 7.1. For nondegenerate weak detonations, D(0) = 0, while 

(7.3) D'(0) = (dd/ds) (u+, s) |a_| < 0, 

where d(u+,s) is the Melnikov separation function defined in (3.5), dd/ds is as 
given in (3.9), and a as usual denotes a — s = df(u) — s. 

Proof. Following the steps of the previous case, we obtain in place of (6.12): 

(u x ^2 y 
|«_| [u 
z x z% 

y,a as defined previously. (Note: the only difference is in the new 2,2 entry la- 
coming from the righthand side of (7.2).) This in turn gives 

(Ux 

z x 

= |a_ I det 

where y := y - (^)w^, o := a - (j^jz^ . But, I = det f U _ x £j , just as 

in the standard, undercompressive shock case treated in [GZ], can be recognized as 
the Melnikov integral (3.9), by expressing I = I~ — I + , where 

u x y± 

Zx ? ± 

with (y + ,a+) = (y+, a+) ^and^(y", a~) = (y",cr _ ) - ([it]/|a_|)(u2 and ob- 

serving that I ± (since do y^, a ± ) satisfy the same ODE 

l'-(a + k<p)I = det(^ (u - u + ) + q(z - z + )\ 

\z x (z x /s) J u=q 

and using Duhamel's principle to express I^ 1 as the integral from ±00 to of the 
integrand in the first line of (3.9); see, e.g., [GZ] for details. The evaluation in the 
second line follows from the identity — s(u — u+) — sq(z — z+) = v! — (f(u) — f(u+)), 
a rearrangement of the traveling-wave ODE, and elementary calculation. The sign 
of the integral then follows from z x > 0, u x < 0, and /' > 0. B 



u 2 




V 


\a- 







z 3 




a 


Ux 


y 


), 


Zx 


a 
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Corollary 7.2. For nondegenerate weak profiles, there holds 



(7.6) 



T= sgn \a-\(dd/d(u+, s))u x (+oo)z x (-oo) > 0, 



consistent with stability. 

Proof. The result of Proposition 6.2 is independent of the case, and so the 
evaluation (6.13) of sgn D(+oo) remains valid, as does the homotopy to A = 



0. However, det \ - _ ) now becomes asymptotic to the explicitly prescribed 



giving a combined result of sgn D(+oo) = sgn u x (+oo)z x (— oo) < 0. Combining 
with our previous calculation of sgn D'(0) < 0, we obtain the result. B 

Remark 7.3. As in [GZ], notice the relation between the key Melnikov integral 
dd/ds and well-posedness of the Riemann problem, Remark 4.2. 



Our methods, in the basic form presented here, require spatially exponential 
decay of the profile, so do not apply to the remaining class of degenerate, subalge- 
braically decaying profiles. 4 Nonetheless, we can make one or two comments of a 
general nature regarding their stability. 

Degenerate detonation profiles. Degenerate detonation profiles occur only in 
the special case u+ = Ui, and generically do not appear in Riemann solutions; in 
this sense, they may be considered a technical curiosity arising from the physically 
artificial cutoff at Ui of the ignition function (for further discussion, see Section 9). 
When they do appear, they occur as a one-parameter family of degenerate strong 
detonation profiles, bounded above by a nondegenerate, and presumably stable, 
strong detonation profile, and below by a nondegenerate gas-dynamical shock profile 
followed by a degenerate weak detonation profile. 

This is somewhat reminiscent of the one-parameter families of overcompressive 
shocks studied for nonstrictly hyperbolic conservation laws in [Fre.2,L,FreL]. As in 
that situation, degenerate strong detonations cannot be singly orbitally stable, but 
only as a family; moreover, in the small viscosity limit (i.e., the true ZND limit, 
without rescaling), there is the same situation that the L 1 difference between profiles 
corresponding to fixed orbits goes to zero, so that stability if it exists at all cannot 
be expected to be uniform with respect to L 1 in the small viscosity limit. (Note 
that, since subalgebraically decaying, they are infinitely far in L 1 from the bounding, 
exponentially decaying strong detonation.) On the other hand, conservation of mass 




endstates 



det 




Section 8. Stability of degenerate waves. 



4 See however [H.1-2,HZ,SS] for extensions to more general situations. 
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(in this case, of f(u + qz)dx) no longer determines uniquely the time-asymptotic 
profile, and it is not clear that any degenerate profile should be expected to be 
stable under perturbation even for fixed viscosity. 

In any case, at the level of Riemann solutions, degenerate strong detonations 
cannot be distinguished from the bounding strong detonation, so give no contradic- 
tion with the conclusions of [CF] based on the ZND model. On the other hand, as 
we saw in Section 4, existence of degenerate weak detonation profiles does lead to 
nonuniqueness in solutions of the Riemann problem for the special case u+ = U{. 
One might hope, therefore, that they could be discarded on grounds of instability, 
thereby validating ZND conclusions for k sufficiently small. 

Degenerate deflagration profiles. Contrary to the situation of degenerate deto- 
nations, weak deflagrations are always of degenerate type, and play an important 
role in the solution of Riemann problems. In particular, stability of (degenerate) 
weak deflagration profiles is necessary for any k > in order that there exist a 
general Riemann solution composed of stable viscous profiles, hence expected from 
the ZND point of view. Thus, stability vs. instability of these waves has important 
philosophical consequences, and is deserving of further study. We point out only 
that, discarding cross terms in the linearized equations, all remaining terms for 
weak deflagrations are favorable for a basic L 2 energy estimate, in contrast to the 
detonation case: the term (au) x , because a x > (expansivity) due to monotonicity 
of the wave, and the term qkdipzu because dtp < in the deflagration regime. 

The stability of degenerate waves of either type (detonation or deflagration) is 
likely to be sensitive due to the subexponential rate of decay of the background 
profile, as for example in the analogous case of KPP waves; see, e.g., [He]. 



Finally, we discuss various elaborations that can be accomodated in the basic 
model without affecting the analysis: 

Multi-species and multiple reactions. Combustion involving more than one 
reactant or reaction can be modeled abstractly, following [FD], by the use of progress 
variables 



each < Aj < 1 denoting the progress toward completion of a single idealized reac- 
tion. As described in [FD], the same framework can be used to describe arbitrarily 
complex reactions in a compact form,with the Xj now representing linear combi- 
nations of progress variables for several simpler reactions. The progress variables 
satisfy a rate equation 



modeling the composite chemical reactions, and enter the kinematic equations via 



Section 9. Extensions. 



(9.1) 



A — (Ai, A2, • • • , A m ), 





(9.3) 



(u - QA) t + (f(u) = u xx , 
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where Q = (q±, q 2 , • • • q m ) denote heats of reaction, qj 
for qj < for endothermic reactions. 

For comparison, consider first a single reaction A 
the present notation by a single progress variable Ai, 

(9.4) \ 1 = k(\ 1 ,u)(l-\ 1 ), 
a typical physical rate k being the Arrhenius rate 

(9.5) k(X-i_,u) = e-^, 

where T is temperature, R is the gas constant, and Eq is activation energy: the 
threshold determining ignition temperature. In the simple model considered previ- 
ously, z corresponds to (1 — Ai), and kip to the rate function k(-, •). Modification 
of the Arrhenius rate by a low-temperature cutoff Tj is a standard device in the 
combustion literature, circumventing the "cold boundary" problem that the Arrhe- 
nius kinetics do not permit a stationary unburned state and therefore preclude the 
existence of traveling wave profiles. 

A more realistic model of combustion consisting of two, consecutive reactions, 
A — > £?, B — > C, is modeled by the rate equations 

(96) r X 1 = k 1 (X,u)(l-X 1 ), 

I A 2 = k 2 (X, u)(Ai - A 2 ), 

with either (i) qi, q 2 > (exothermic/exothermic), or (ii) q\ > 0, q 2 < (exother- 
mic/endothermic). As above, a standard choice of rate function would be kj = 
ipj(u)kj, with kj constant, and <pj ignition functions of the usual form. As pointed 
out in [FD], the traveling wave structure in case (i) is very similar to that for a 
single reaction; on the other hand, (ii) introduces important new phenomena at 
the ZND level, explaining e.g. an observed experimental shift from CJ to weak 
detonation in ignition problem. See [FD], p. 168-173, for further discussion. Three 
or more reactions can easily be included, but do not appear to result in significant 
new phenomena [FD]. 

It is straightforward to show that the main conclusions of the previous sections 
remain valid for equations (9.6), or indeed any system of reactions having such a 
triangular structure: more generally, for any reaction system in which A = (1, • • • 1) 
is stable above the ignition temperature (s) as a rest point of rate equation (9.6). In 
particular, we recover the result for strong detonations that the stability index is al- 
ways positive, consistent with stability. Likewise, we recover the same expression for 
r in the case of weak detonations, but now the Melnikov integral for (dd/ds)( u s ) is 
more complicated and its sign (apparently) no longer explicitly evaluable. It would 
be very interesting to investigate structure and stability of traveling waves for the 
2-reaction model above, especially in light of anomalous behavior (e.g., a shift in 



> for exothermic reactions, 

— > B. This is represented in 
with rate equation 
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the ignition problem from CJ to weak detonation) predicted in [FD,B] in the the 
ZND setting, in the high-activation energy limit. 

Reaction-dependent equation of state. As pointed out in [CHT], a more 
realistic assumption is that the gas-dynamical equation of state (EOS) / depend 
not only on u, but also on z, through the chemical makeup of the gas. With this 
change, the linearized equations (5.1), (5.2) become 

^ ^ ( u t - q{kip'(u)uz - k<p(u)z) + (au) x + (f3z) x = u xx 

\ z t — sz x = — kip'(u)uz — kip(u)z. 



^ ^ ( u t + qz t + (au) x + (pz) x - sqz x = u xx 

\ z t — sz x = — k(f'(u)uz — k(f(u)z, 

where a := f u (u, z)(x) — s and (3 := f v (u, z)(x). It is straightforward to check that 
all of our stability analysis goes through as before in this more general setting, to 
yield exactly the same geometric stability conditions. In particular, T > holds 
always for a strong detonation profile if it exists. Likewise, the expression derived 
for r in the weak detonation case, including our calculation of the Melnikov integral 
dd/ds, remains valid, though the sign of the integral may change depending on the 
details of the equation of state. 

Regarding the existence problem, and the related question of stability of weak 
detonation profiles, we are led naturally to a simple condition, 

(9.9) /, < 0, 

under which all of the conclusions of the paper generalize to the case of a reaction- 
dependent EOS. For, this implies that F z = f z — qs < 0, where F = as in (3.4) 
defines the nullcline u' = for the traveling-wave equations, and so the nullcline is 
again a graph over z. This was essentially the only property used in the qualitative 
analysis of the phase portrait, and so we recover all related conclusions, including 
the important property of monotonicity in u and z of weak detonation profiles. 
Likewise, (9.9) implies that f(u, z) — f(u + , z+) > along weak detonation profiles, 
since f u > 0, f z < and u and z are, respectively, monotone decreasing and 
increasing. Thus, we obtain the key monotonicity in s asserted in Proposition 3.3, 
now in the restricted case z < 1, in particular, the conclusion that dd/ds < at 
a profile d = 0, along with the other stated monotonicity results (the arguments 
for which are unaffected by dependence of / on z), and thereby the remaining 
conclusions of the paper. 

On the other hand, when the monotonicity condition (9.9) is violated, we see 
no reason why anomalies in both existence and stability theory might not occur, 
and this seems an important direction for further investigation. In the full, reacting 
Navier-Stokes equations, the flux functions, taking gas composition into account, 
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are just the convex averages, weighted by z, of corresponding fluxes for pure un- 
burned and pure burned gas, and similarly for the coefficient c v relating temperature 
to internal energy. We conjecture that the corresponding monotonicity condition 
in this physical setting reduces to the condition that the gas constant 7(1) of the 
unburned gas be less than the gas constant 7(0) of the burned gas, assuming that 
each separately obeys the EOS of an ideal, polytropic gas: i.e., from the kinetic 
theory of gases point of view, the average number of internal degrees of freedom 
per molecule n decrease upon reaction (recall: 7 = (n + 2)/n; see, e.g., [Ba], pp. 
37-45). The derivation of an analog of condition (9.9) by analysis of the ZND limit 
we regard as another very interesting open problem. 
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